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ABSTRACT 

Numerical approximate computation can solve large and 
complex problems fast. It has the advantage of high effi- 
ciency. However it only gives approximate results, whereas 
we need exact results in many fields. There is a gap be- 
tween approximate computation and exact results. A bridge 
overriding the gap was built by Zhang, in which an ex- 
act rational number is recovered from its approximation 
by continued fraction method when the error is less than 
1/((2JV + 2)(iV - 1)N), where N is a bound on absolute 
value of denominator of the rational number. In this pa- 
per, an improved algorithm is presented by which a exact 
rational number is recovered when the error is less than 
l/(4(iV - 1)N). 
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1. INTRODUCTION 

Numerical approximate computations have the advantage 
of being fast, flexible in accuracy and being applicable to 
large scale problems. They only give approximate results 
and are applied in many fields. However, some fields such 
as theorem proving, need exact results and symbolic com- 
putations are used to obtain the exact results. Symbolic 
computations are principally exact and stable. They have a 
high complexity. They are slow and in practice, are appli- 
cable only to small systems. In recent two decades, numer- 
ical methods are applied in the field of symbolic computa- 
tions. In 1985, Kaltofen presented an algorithm for perform- 
ing the absolute irreducible factorization, and suggested to 
perform his algorithm by floating-point numbers, then the 
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factor obtained is an approximate one. After then, numer- 
ical methods have been studied to get approximate factors 
of a polynomial 4 11 15 16 17 18.. In the meantime, nu- 
merical methods are applied to get approximate greatest 
common divisors of approximate polynomials !l 2- [12] [5] , to 
compute functional decompositions^, to test primalitv|10| 
and to find zeroes of a polvnomial[13|. In 2000, Corless 
et al. applied numerical method in implicitization of para- 
metric curves, surfaces and hypersurfaces[3]. The resulting 
implicit equation is still an approximate one. 

There is a gap between approximate computations and ex- 
act results [i"9]. People usually use rational number compu- 
tations to override the gap [7]. In fact, these are not approxi- 
mate computations but big number computations, which are 
also exact computations. In 2005, Zhang et al proposed an 
algorithm to get exact factors of a multivariate polynomial 
by approximate computation [20] but they did not discuss 
how to override the gap. In [9], Cheze and Galligo discussed 
how to obtain an exact absolute polynomial factorization 
from its approximate one, which only involves recovering an 
integer from its approximation. They did not discuss how 
to obtain an exact rational number from its approximation. 
Command convert in maple can obtain an approximate ra- 
tional number from a float if we set variable Digits to a 
positive integer. When variable Digits is taken to different 
positive integers, a different rational numbers are obtained. 
Which one is the rational number we want? We do not know. 
In some cases, what's more, none of the rational numbers is 
the one we want. For example, if we want to get 1/7, we 
take 0.1196013289 as its approximation. However, no mat- 
ter what taking variable Digits to, we can not obtain 1/7 by 
command convert. In [22], Zhang and Feng systematically 
discussed how to obtain the exact result from its approxi- 
mation. They proved that the exact rational number can be 
obtained from its approximation when the error is less than 
1/(2N(N — 1)), but in practice, the algorithm requires that 
the approximate error is less than 1/((2K + 2)N(N — 1)), 
where N is an bound on absolute value of denominator of 
the exact rational number and K > N. In this paper, we 
propose an improved algorithm which can recover the exact 
rational number from its approximation when the error is 
less than 1/(4N(N - 1)). 



The remainder of the paper is organized as follows. Sec- 
tion 2 gives a review of some properties of continued frac- 



tion, which are used to prove our theorems later. Section 
3 proves that the rational number we want is the only one 
satisfying the error control e < 1/(2N(N — 1)), and then 
shows that a kind of rational numbers can not be obtained 
from its approximation by continued fraction method when 
the error l/(4N(N - 1)) < e < 1/(2N(N - 1)); and finally, 
an improved algorithm is proposed. Section 4 gives some 
experimental results. The final section makes conclusions. 



Theorem 3. Each convergent is nearer to the nth con- 
vergent than any of the preceding convergents. In symbols, 
if the nth convergent ts taken to be [ao; a ly a-z, ■ ■ • , a n ] = x, 
then 

|[ao; ,ai,a 2 , ■ ■ ■ ,a r ]—x\ > |[a ; ,ai,a 2 , ■ ■ ■ ,a s ] — x\ 
for all r < s < n. 



2. PROPERTIES OF CONTINUED FRACTION 

A continued fraction representation of a real number is one 
of the forms: 



ao + 



ai + 



(1) 



a 2 + : 



'3+' 



where ao is an integer and 01,02,03, •■ ■ are positive inte- 
gers. One can abbreviate the above continued fraction as 
[oo;ox, 02) • • • ]. For finite continued fractions, note that 



a ; ai, a 2 , 



, a„, 1] = [ao; oi, a 2 , 



, a„ + 1] 



So, for every finite continued fraction, there is another finite 
continued fraction that represents the same number. Every 
finite continued fraction is rational number and every ra- 
tional number can be represented in precisely two different 
ways as a finite continued fraction. The other representation 
is one element shorter, and the final term must be greater 
than 1 unless there is only one element. However, every 
infinite continued fraction is irrational, and every irrational 
number can be represented in precisely one way as an infinite 
continued fraction. An infinite continued fraction represen- 
tation for an irrational number is mainly useful because its 
initial segments provide excellent rational approximations to 
the number. These rational numbers are called the conver- 
gents of the continued fraction. Even-numbered convergents 
are smaller than the original number, while odd-numbered 
ones are bigger. If successive convergents are found, with nu- 
merators hi,h,2, ■ ■ ■ , and denominators fei,fca, • •• , then the 
relevant recursive relation is: 

hn = CLnhn — 1 ~\~ h n — 2, kn = &nkn — l ~\~ kn — 2- 

The successive convergents are given by the formula 

+ h 

n — 2 

k n Q>nk n —1 ~\~ k n — 2 

where h-i — 1, h-2 = 0, fc-i = and k-2 = 1. Here are 
some useful theorems [8]: 

Theorem 1. For any positive x £ R, it holds that 

xhn-i + h n -2 



[ao, ai, 



-i,x = 



xk n —i -j- kn—i 



(2) 



Theorem 2. The convergents of [ao, 01, 02, • • • ] are given 



by 



and 



[a ,ai, • ■ ■ ,a„\ = — 

fcn 



knhn — 1 kn — lhn — ( 1) • 



Proof: Denote x = [ao; 01, 02, • • • ,a r ,x r +i], where av+i 
[a r +i, • • • , On]. From theorem [TJ it holds that 

x r+ ih r + h r -i 



Xr-^-^ky ~\~ kf — 1 

Accordingly, we can deduce as follows: 

X(x r +lk r + fc r _l) = Xr + lhr + h r -\ 

— r" Xr-^-\{^xky — [xky — ]_ hy — 1) 

Dividing above equation by x r +ik r yields 

h r , k r —\ . , h r —\ . 

x ~ t = — r^ x ~ t — ) 

Since x r +i > 1 and k r > k r -i > 0, we have 
0<(^)<1 

X r -\-\K, r 

Therefore, it is proved that 

1 hr 1 1 hr — 1 I 

I X -F I<|X ~-W '■ 

The proof is finished 

Theorem 4. Let x = [ao, 01,02, • • •] and h n = a n h n -i + 



and 



n kn — 


1 + kn-2- Then it holds that 




hn h n —± 


1 




kn kn—l 


kn kn—l 



kn(kn-\-l ~\~ kn) 



hn 
kn 



kn kn-\ 



3. AN IMPROVED ALGORITHM FOR RE- 
COVERING THE EXACT NUMBER FROM 
ITS APPROXIMATION 

In this section, we will solve such a problem: Someone has a 
positive rational number ^ in his mind, and you only know 
an upper bound N of denominator of the rational number, 
and he can not tell you the rational number but an approxi- 
mation of the rational number at any accuracy. How do you 
compute the rational number from one of these approxima- 
tions? Let's attack this problem. Without loss of generality, 
we always assume that m,n are positive numbers. At first, 
we have a lemma as follows: 



Lemma 1. m,n,p,q are integer, andpn > 0. If\^ — p \ < 
-2-, then 2H = 1. 



Proof. 



— p | = p„ ■ Noticing |pm — qn\ is a non- 
negative integer, and |22 — 2,| < yields \pm — gn| < 1. 
Hence \pm — qn\ = 0. That is 2 = £. The proof is finished. 



So, we have r = [0, 1, n — 2, 2, n, n — 1] . It is obvious that we 
can not obtain n—l/n from [0, 1, n — 2, 2, n, n— 1]. Therefore, 
we need smaller neighborhood so as to recover the exact 
rational number by continued fraction method. 



Corollary 1. m,n,p,q are integers and p > 0,n > 0. 
ZrfiV>max{p,n,2}. If\^-l\ <w ^, then 2 = f. 



Proof: When p 7^ n, it holds that pn < AT(iV — 1). Hence 
I — — - 1 < ... * . , < — . According to lemma [0 it is 

1 n p N(N — 1) — pn ° 1 — • 

obtained that — = -. When p — n, we have 

n p r 

m g ! |m — q\ I 



n p 1 n N(N-l) 



\m — q\ < 



N(N — 1) ~ N -1 



< 1 



So, it holds that m — q. The proof of the corollary is fin- 
ished. 



And now, we discuss how to obtain rational number from 
its approximation by continued fraction method. Let ni /n\ 
be a rational number and ro its approximation. Their con- 
tinued fraction representations are 712/ ni — [ao,<Ji,-' - , Oi] 
and ro = [60, fei, • • • , bjvr] respectively. We wish that a,i = bi 
for £ = 0, 1, 2, • ■ ■ , L — 1 and for the last term of the con- 
tinued fraction representations of n,2/ 'ni, either ol = 6l or 
sl - 1 = 6_l, so that we can get nijn\ from [60, 61, • • • , 0z,+i]. 
This is the following theorem: 

Theorem 6. Let 112/ n\ be a rational number and ro its 
approximation. Assume that n2 ,n\ are coprime positive num- 
bers, where n2 < n\,and m > 1. The representations of 
n2/ni and ro are [ao, fli, • • • ,ol] and [bo, 61, • • • , 6a/] respec- 
tively. If |ro — W2/ni| < l/(4rti(rti — 1)), then one of the 
following statements must hold. 



Theorem 5. Letx= — be a reduced proper fraction, and 
N > max{n, 2}. Assume'that \x - w\ < l/(2N(N - 1)). // 
we get positive rational number p/q such that \p/q — w\ < 
1/(2N(N - 1)) , where q < N , then it holds that x = q/p. 



• ai = bi (i = 0, 1, • ■ • , L); 

• at = bi (i = 0, 1, ■ ■ ■ , L — 1), ol — 1 = fee, and bL+i — 
1. 



Proof: From the assumption of the theorem, we have \x — 
q/p\ < l/(N(N-l)). According to corollaryQ] it holds that 
q/p = m/n = x. The proof of the theorem is finished. 

Theorem [5] shows us as follows. One has a rational number 
— in his mind, if he tells you an approximation w such 
that — - w\ < 1/(2N(N - 1)), then there is an unique 
rational number whose denominator is less than iV in the 
neighborhood (w - 1/{2N(N - l)),w + 1/(2N(N - 1))). 

The remaining question is as follows. How do we fetch out 
the unique rational number in the neighborhood? We wish 
get it by continued fraction. Unfortunately, we can not al- 
ways fetch out the unique rational number in the neigh- 
borhood (w-l/(2N(N -l)),w + l/(2N(N -1))) by con- 
tinued fraction method. One counterexample is the ratio- 
nal number such as (n — l)/n for n > 1. Let us show 
this: set N — n, and its approximation r := (2n + 2n 3 — 
4n 2 — l)/(2n 2 — 2n + l)/n. One can check its error d = 
\/{2n{n - 1) + 1) < 1/(2N(N - 1)). However, one can not 
recover rational number (rt — l)/rt from its approximation 
by continued fraction method. In fact. First we can easily 
compute (n — l)/n = [0, l,n — 1], continued fraction repre- 
sentation of n — l)/n. And then compute continued fraction 
representation of r as follows. 



2n + 2n 6 - 4n 2 



1 



1 



2n 2 



n+1 



ri 



i-2 



r 3 



(2n 2 - 2n + l)ra ~" r 1 + 2n + 2n 3 - 4n 2 - 1 

2n 2 - n + 1 
2n + 2n 3 - An 2 - 



1 1 + n z - n 
— n _ 2 + I 

ri 2n 2 —n + 1 



1 + n z 



2n 2 - n + 1 

71-1 1 



1=2+ 

r2 1 + n 2 — n 



1 



According to assumption of ri2 < m, we have that ao = 
0, and bo = 0. Hence ao = feo- In order to finish the 
proof of theorem [(J we need two lemmas. Due to n%jn\ — 
[ao,ai, ••• , cll] and ro = [60, &i, • • • ,6a/], we have the fol- 
lowing expansions: 

m 



and 



n 3 n 2 , n 4 
— = aiH , — = a 2 -\ , 

n,2 ri2 713 ri3 

riL-i 1 

= OL-l H , TIL = a I, 

TIL TIL 



1 1 

— = 61 + n, — = 02 + r 2 , 
ro ri 

1 . 1 

: ol + r L , 



(3) 



6m (4) 

r_L-i r A /-i 
Denoting di = ri~ ni+2/n.i+i, we have a lemma as follows: 

Lemma 2. iet 712/ni 6e a rational number and ro its ap- 
proximation. Assume that ri2,n\ are coprime positive in- 
tegers, where 712 < n\,and n\ > 1. T/ie representations 
ofriz/ni and ro are [ao, ai, ■ ■ ■ , az,] and [feo, 61, • • • , &m] re- 
spectively. And assume that ai = 6i /or i < k < L(k is a 
positive inteqer). Then when Idtl < rr-, it holds 



f/ia£ a fe+ i = b k+1 for k < L - 1; when |dz,_i| < „ L(nL+1) 
holds that ol — bL or aL — 1 = 6l • 



. it 



Proof: At first, we show that under the assumption of the 
lemma if we have 

nl + id k 



nu+2{nk+2 + n k+1 dk) 



< 



llk+2 



(5) 



1 + n 2 



- n 



r?, 



then, it holds that ak+i = bk+i for k < L — 1, and ak+i = 
bk+i or ak+i — 1 = 6fc + i for k = L — 1. We discuss it in two 



Case l(k < L — 1): From d k = r k — rik+2/nk+i, it holds that 
i"k — nk+2/n.k+i + d k . Hence we have that 



1 

rk 
1 



7ifc+i 
rik+2 
nfc+i 



n k+2 {n k+ 2 +n k+1 d k ) 
n 2 k+1 d k 



r k n k+2 n k+2 (n k+2 + n k+ id k ) 

=>■ — 1fe+l H 7 j TT 

r k n k+2 n k+2 (n k+ 2 + n k+1 d k ) 

= b k+1 + r k+1 

Hence, obviously, a k +i = b k+ i if and only if 

n 2 k+1 d k 



< 



n k +3 

n k+2 n k+2 (n k+ 2 + n k+1 d k ) 



< 1. 



(6) 



Therefore, if inequality ([5]) holds, then above inequality is 
guaranteed. 

Case 2: (when k = L — 1) We have 



1 



= a L 



n 2 L d L -! 



n L +i(n L +i + n L d L -i\ 
From the above equation, if 



= b L + r L 



n h dh- 



< 



n i+ i(n i+ i +n L rf L _i) n L+2 

then aL = &l for dz,_i < 0, and ol — 1 = &l for dL-i > 0. 
Therefore, we have shown that if inequality ((5| holds, then 
Ofc+i = b k +i for k < L — 1, and either afc+i = b k +i or 
afc+i — 1 = b k+1 for fc = L — 1. 



On the other hand, we have 



t+i|dfc| 



\n k+2 (n k+2 + n k+ id k )\ 
n k +i\d k \ 



n k+ 2{n k+2 + n k+ id k ) 
_ n k+ i\d k 

- n k+2 (\^+d k \) ~ n k+2 \(^ - \d k \)\ 

So, in order to ensure inequality JSJl, we only need it holds 
that 



n k+1 \d k \ 



n k+ 2\(^-\d k \)\ 
Solving inequality @ yields 

n k+ 2 



< 



n k+ 2 



\d k \ < 



(7) 



(8) 



n k+ i(n k+1 + 1) 

When k < L — 1, we have that n k+ 2 > 1. So, it holds that 

1 < n k+2 

n k+1 (n k+ i — 1) ~ n k+1 (n k+1 + 1) 

Accordingly, it is obtained that 

1 



\d k \ < 



n k +i{n k -\ 



1) 



(9) 



When k = L — 1, we have that ul+i = 1, so it is obtained 
that 



\d L -i\ < 



nL{n L + 1) 



(10) 



Lemma 3. Let 712/721 be a rational number and ro its ap- 
proximation, where n 2 ,n\ are coprime positive integers, and 
712 < Tii, and ni > 1 . The continued fraction representations 
0/7x2/711 and ro are [ao, ffli, • • ■ , ol] and [60, 61, •• • ,6m] 7~e- 
spectively. Denote di — — 7i i+ 2/7i i+ i /or i = 0, ■ ■ ■ , L. 
Assume that a; = 6^ /or i < fc < L — l(k is a positive integer 
). Then when \d k \ < t-t, it holds that 



I , n k+1 {n k+1 - 1) 

"fc+i < 7 7T °* 

n k+ 2\n k+ 2 — 1) 



(11) 



Proof: Under the assumption that o^ = fe; for i = 0, 1, • • ■ , fc, 

2 j 

from equation ©, we get dfc+i — fc+1 k 

we deduce a relation as follows: 



Idfc+il = 



n k+ i\d k \ 



n k + l{ n k + 2+ n k+l d k) 

n k+ i\d k \ 

n l+2 + n k+1 n k+2 d k nfc+2 (!^±l + d k ) 

n k+ i(n k+1 — l)\d k \ 

«w( ! ^ + K + i-i)i) 



Hence 



71fc+l(71fci 



D\d k 



n k+2 {n k+2 - 1 + ""^ 7° +2 + (n-k+i ~ 1)4) 



When Idfel < 



. it holds that 



"fc + l~"fc + 2 



+(n fc+ i- 



l)d k > 0. Hence we have a relation between dfc+i and d k : 
I n fc+ i(n fc+ i - 1) 

fflfc+l < 7 rv Ofc 

n k+ 2(n k+ 2 — 1) 
The proof of the lemma is finished. 

And now, let us prove the theorem. If \do\ = |fo — n2/ni\ < 
l/(4ni(m — 1)), From lemma[3] we can get 

\dt\ < l/(4n j+ i(n i+ i - 1)) 

for i = 0, • • • , L — 1. Note that ul > til+i = 1 and 

1 1 < 1 

4n i+ i(n i+ i - 1) n i+ i(n i+ i + 1) n i+1 (n i+ i - 1) 

when m+i > 1. So, it holds that 

1 1 



< 



4n i+ i(n i+ i - 1) n i+1 (n i+1 + 1) 



for i = 0, • ■ ■ , L — 1. According to lemma[2] the proof of the 
theorem is finished. 

For an unknown rational number 712/711 and its approxi- 
mation ro, theorem [S] shows that 712/ni = [60, • ■ • ,6^,] or 
712/711 = [bo, 61, ■ ■ ■ , b L , 1] when |r - n 2 /ni| < l/(4ni(rai - 
1)). However, for practical purpose, we hope the restriction 
on ni > 1 and m > 712 can be lifted. So we have following 
theorem: 



Theorem 7. Letno/ni be a reduced rational number and 
r its approximation. Assume that 7io,7ii are positive integers 
and N > max{7ii,2}. The continued fraction representa- 
tions o/tio/tii and r are [ao, 01, •• ■ ,cll] and [bo, 61, ■ ■ ■ , 6m] 
respectively. If\d\ = |j — no/ni| < 1/(4^(^ — 1)), then one 
of the following two statements must hold 



The proof of lemma [2] is finished. 



• Oi = bi for i = 0, • • • ,L; 



• a,i — bi for i = 0, • • ■ , L — 1, and ol = fli — 1, = 
1. 



Proof: We prove the theorem in three cases: 

Case 1 (m > 1, n < m): From 1/(4A(7V-1)) < l/(4ni(m- 

1)) and theorem [6] the theorem holds. 

Case 2(ni = 1): We have that ao = rio/ni. If d = r — 

no/ni > 0, then bo = ao and ro = r — ao < 1/(4 x 2(2 — 1)). 

If d = r — no/ni < 0, then 

r = ao — \d\ = ao — 1 + 1 — \d\ = bo + 1 — |d| 
=>■ ro = r — bo = 1 — |d| 



On the other hand, we have 

|d| < 1/8 => 2|d| < 1 ^ |d| < 1- |d| 

^ < -^L < 1 6r = 1 
1 — l"l 

So, we have that f>o = ao — 1, &i = 1. 

Case 3(no > ni > 1): From no/m = ao + n^/m, it holds 
that no/m — ao = n%/m. On the other hand, we have that 
|n /ni — r| < 1/(4N(N - 1)) < l/n%. So, we can deduce 
that a < r < a + 1. Accordingly, it holds that & = «o- 
Hence, we have 

|d| = |r — n /ni| = |6 +r - a - n 2 /m = \r - n 2 /m 
= do < 1/(4A(A-1)) < l/(4ni(m - 1)) 

And now, we have m > 1 and 712 < m, which is case 1. 
Therefore, the proof is finished. 

For simplicity, set s — L when [ao, ai, • • • , ax] = [fco, 61, • • • , &t] 
and s = L + 1 when [ao, 01, ■ • • , a L ] = [b , 61, ■ ■ ■ ,b L , 1]. 
For an unknown rational number no/ni and its approxi- 
mation r, theorem [7] shows that no/m = [bo, ■ • ■ ,b s ] when 
\r — no/ni| < 1/(4N(N — 1)). However, we do not know 
what the number s is. The following theorem shows us how 
to get [6 , ■ ■ ■ , b 3 ]. 



Theorem 8. Letno/ni be a reduced rational number and 
r its approximation. Assume that no,ni are positive integers 
and N > max{m, 2}. We have r = [60, bi, • • • , bit] and \r — 
n /ni| < l/(4N(N - 1)). Denote no/m = [&o,&i,--' ,&»]- 
Then, for any positive integer s < t < M , the denominator 
of rational number g — [bo, bi, • ■ ■ ,bt] is greater than N . 



Proof: Proof is given by contradiction. Denote by m the 
denominator of g. Assume m < N. From theorem[3] it holds 
that [r - [60, 61, •• ■ , 6*]| < |r — [60, 61, •• • , 6 S ]|. Noting that 
|r- [60, ,MI = \r-n /m\ < 1/(4N(N - 1)) yields 
\r- [bo,h, ■ ■ ■ ,b t ]\ < 1/(4N(N- 1)). According to theorem 
[5] it should hold that [60, 61 , • • • , b s ] = [60 , 61 , • • • , bt] ■ This 
contradict to that t > s. The proof is finished. 

Based on theorem[7]and theorem[5l an algorithm for obtain- 
ing the exact number is as follows: 



Algorithm 1. Input: a nonnegative floating-point num- 
ber r and a positive number N ; 
Output: a rational number b. 

Step 1: Set i = 0, tem = r, h-\ = 1, h-2 = 0, 
k-i = 0, and k-2 ~ 1; 

Step 2: Get integral part of tem and assigning it 
to a, assigning its remains to b. 

Step 3: Compute hi = a * hi-i + hi-2 and ki = 
a * ki-i + ki-i. If h > N,then goto Step\Bj 

Step 4: Set i := i + 1; 

Step 5: Set tem = | and goto Step 0- 

Step 6: Computing hi-i/ki-i and assigning it to 
b. 

Step 7: return b. 

The correctness of algorithm [T] is obvious from theorem [7] 
and theorem [5] 

4. EXPERIMENTAL RESULTS 

The following examples run in the platform of Maple 10 and 
PIV 3.0G, 512M RAM. They take little time for obtaining 
exact rational numbers from their approximations, so we do 
not show time. 

Example 1. Let a be unknown rational number. We only 
know a bound of its denominator N — 170. According to 
theorem [7J Computing rational number a as follows: Com- 
pute d = 1/(4 * JV * (jV - 1)) = 1/114920. Assume that 
we use some numerical method to get an approximation 
b = .8106421859 such that \a - b\ < 1/d. Calling algorithm 
Q] yields as follows. 

4 13 17 30 107 137 518 
' ' 5' 16' 21' 37' 132' 169' 639 
When algorithm[JJfinds that the denominator of ||| is greater 
than N, it outputs 

Example 2. Let a be unknown rational number. We only 
know a bound of its denominator N = 1790. According to 
theorem [7] Computing rational number a as follows: Com- 
pute d = 1/(4 * N * (N - 1)) = 1/12809241. Assume that 
we use some numerical method to get an approximation 
b = .178870799516605 such that \a - b\ < 1/d. Calling 
algorithm [JJ yields as follows: 

1 1 2 5 17 22 149 171 320 1131 

' 5' 6' IT' 28' 95' 123' 833' 956' 1789' 6323 

When the algorithm finds that the denominator of 1131/6323 
is greater than N, it outputs 320/1789. 

Example 3. Let a be unknown rational number. We only 
know a bound of its denominator JV = 18. According to the- 
orem [7] Computing rational number a as follows: Compute 
d = 1/(4* AT* (N- 1)) = 1/1225. Assume that we use some 
numerical method to get an approximation b — 1.881536615 
such that [a— 6| < 1/d. Calling algorithm [T] yields as follows. 

1,2,15/8,32/17,111/59 



When the algorithm finds that the denominator of 111/59 
is greater than N, 32/17. 

Example 4. This example is an application in obtaining 
exact factors from their approximations. Let p — — 16 — 
56y - 48z + 64x 2 - 32xy + 4Sxz - 45y 2 - 96yz - 21 z 2 be a 
polynomial. We want to use approximate method to get its 
exact factors over rational number field. First, we transform 
p to a monic polynomial as follows: 

2 1 3 45 2 3 27 2 7 3 1 

p = x xy H — xz y yz z y z 

' 2 y 4 64 tf 2 y 64 8 y 4 4 

the least common multiple of denominators of coefficients of 
polynomial p(x,y,z) is 64, which is an upper bound [21] of 
denominators of coefficients of the monic factors of polyno- 
mial p. Taking N = 65 yields d = 1/(4 * 65 * 64) = 1/16128. 
We use numerical methods to get its approximate factors as 
follows [21]: 

gi = 1.0000a;+.6250000000067y+1.124999999530z+.50000 

g 2 = 1.0000a: -1.1250000000152/-. 3749999995480z -.50000 

the error of coefficients of gi and gi is less than d by the 
numerical methods. According to theorem [7] taking TV in 
algorithm [T] we obtain two exact factors: 

5 9 1 
91 =x+-y+-z+- 

9 3 1 

02 = x y z 

y 8 y 8 2 

5. CONCLUSION 

In [22], a bridge overriding the gap between approximate 
computation and exact results was built. However, the al- 
gorithm in [22] requires the error between approximation 
and exact result is less than l/((2iV + 2)N(N - 1)). In this 
paper, we propose an algorithm that only requires the error 
is less than 1/(4N(N — 1)), which decreases the cost of com- 
puting approximation. Just like the algorithm in [22], our 
method can be applied in many aspect, such as proving in- 
equality statements and equality statements, and computing 
resultants, etc. Thus we can take fully advantage of approx- 
imate methods to solve larger scale symbolic computation 
problems. 
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